Doped Mott insulator as the origin of heavy Fermion behavior in LiV204 
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We investigate the electronic structure of LiV204, for which heavy fermion behavior has been ob- 
served in various experiments, by the combination of the local density approximation and dynamical 
mean field theory. To obtain results at zero temperature, we employ the projective quantum Monte 
Carlo method as an impurity solver. Our results show that the strongly correlated ai g band is a 
lightly doped Mott insulator which -at low temperatures- shows a sharp (heavy) quasiparticle peak 
just above the Fermi level, which is consistent with recent photoemission experiment by Shimoya- 
mada et al. [Phys. Rev. Lett. 96 026403 (2006)]. 
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The discovery of heavy fermion (HF) behavior in the 
3d material LiV 2 04 [l[ was a big surprise since this phe- 
nomenon was previously a hallmark of certain / elec- 
tron compounds. But below a characteristic temperature 
Tk ~ 28 K, the linear specific heat coefficient (7) [![, 
the magnetic susceptibility [l[,the Griineisen parameter 
0, and the quadratic resistivity coefficient Q are also 
for LiV2C>4 extraordinarily large, similar to / electron 
HF compounds and much larger than in other transition 
metal oxides. From 7 ~ 420 mJ/molK 2 , an effective mass 
enhancement of m* /m ~ 25 was inferred Nl. Neutron 
scattering [6|, nuclear magnetic resonance [5| and elec- 
tron spin resonance [7[, as well as muon spin relaxation 
experiments [§] indicate the existence of local magnetic 
moments, which is consistent with a Curie- Weiss suscep- 
tibility |l| in the temperature range 50 to 1000 K. Down 
to the lowest temperatures measured LIV2O4 remains a 
cubic spinel and no long-range magnetic, spin glass, or 
superconducting order was observed. More recently, also 
a sharp Kondoesque peak of width 10 meV was observed 
in photoemission experiments [9( just 4 meV above the 
Fermi energy Ep, with a strong temperature dependence 
similar to that of other HF compounds. This finding is 
supported by the measurement of the magnetic curves 
at T = 1.3 K [l3| which also suggests the existence of a 
sharp peak slightly above Ep. 

The explanation of the HF behavior in LiV204 has 
been a challenge since its discovery. Local density 
approximation (LDA) calculations 0, U, EH show a 



twofold-degenerate and 2eV-wide e g and a nondegen- 
erate a± g band of width 1 eV cross the Fermi energy, 
filled altogether with 1.5 electrons per V ion. This 
LDA bandstructure led one of the present authors (VIA) 
to the proposal 11, [I3J that the a\ g electrons play 
the role of the localized / electrons in conventional 
HF compounds and the e g electrons that of the itin- 
erant valence electrons. On the other hand, the im- 
portance of geometrical frustration originating in the 



spinel structure has been stressed by various authors 
3, [H, EH, E3, EE O 0, ED, H2] , all suggesting different 
explanations for the mass enhancement of LiV 2 04. Nat- 
urally the geometrical frustration suppresses any kind of 
long range order, so that local spin or orbital fluctuations 
should be dominant as suggested in Ref. 20l. 21], 

In this situation, we might expect dynamical mean 
field theory (DMFT) [23} to be good approximation for 
studying the electronic correlations in this material. Re- 
alistic LDA+DMFT [jjcalculations for LiV 2 4 have 
been carried out before [25j , but neglected the a\ g -e g hy- 
bridization, which should be the driving force for the HF 
behavior in the Kondo scenario llj and were further- 
more restricted to temperatures T > 750 K, far above 
Tk- Not surprisingly, a quasiparticle resonance was not 
found and the competition between antiferromagnetic 
direct exchange from the a\ g -a\ g hybridization, ferro- 
magnetic double exchange from the s g s g hybridization, 
and the Kondo effect from the (neglected) aio-e g hy- 
bridization left this LDA+DMFT study inconclu- 
sive. Since that time the more sophisticated projection 
onto Wannier functions has been developed [26( which 
properly takes the orbital off-diagonal hybridization into 
account. Also the problem that conventional quantum 
Monte Carlo (QMC) simulations 27] of the DMFT impu- 
rity problem were restricted to rather high temperatures 
because the numerical effort is proportional to 1/T 3 has 
been overcome by the projective QMC (PQMC) method 
[Si HI E3] forT = 0. 

With these improvements, we reinvestigate LiV204 by 
LDA+DMFT(PQMC) and solve the puzzle why this ma- 
terial shows HF behavior with a sharp Kondoesque res- 
onance above the Fermi level. 

Method. The unit cell of LiV204 contains four V atoms 
and each V atom has three t2 g orbitals, which are split 
into the a\ g orbital and two degenerate e g orbitals due 
to the trigonal splitting. First, we do a LDA calcula- 
tion for LiV204 using the linearized muffin tin orbital 
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FIG. 1: (Color online) Left panel: Band dispersion of the 
effective 2-orbital model (solid line) and total band structure 
(dashed line) of L1V2O4. Right panel: partial ai g and DOS 
for the model. Ef is set to zero. 



basis set [3l|. From this we further construct an effec- 
tive 12 by 12 Hamiltonian by the projection onto Wan- 
nier functions [26]. Since the e g orbitals are degener- 
ate, it is possible to derive a 2-orbital model with an 8 
by 8 Hamiltonian by taking only one of two e g orbitals 
into account. This drastically decreases the computa- 
tional efforts of the LDA+DMFT calculation and hence 
allows for more accurate data. As we will see a posteriori, 
the restriction to one e g orbital will be justified by the 
fact that the e g orbitals play a rather passive role and 
the physics is determined by the a\ g band. The com- 
parison of the band dispersion of this simplified 8-band 
model with total LDA band structure of L1V2O4 in Fig. 
[T] shows that the 2-orbital simplification captures the es- 
sential features of the real compound's band structure. 
It also gives the densities of states (DOS) close to that 
previously reported by LDA, see Fig. 5 in (2o| . 

Second, we supplement this 2-orbital Hamiltonian by 
local intra- ([/) and inter-orbital ([/') Coulomb repulsions 
as well as by Hund's exchange (J), and solve the con- 
structed many-body model by DMFT. It should be noted 
that we explicitly consider the off-diagonal elements be- 
tween e g and a\ g in contrast to all previous calculations 
[25l | . where only the initial LDA aig-e^ hybridization is 
reflected indirectly in the DOS. This a\ g -e^ hybridiza- 
tion is essential for the Kondo effect with localized a ls 
and itinerant electrons 11] . As for the self-energy, we 
only consider the diagonal element, so that the effective 
DMFT impurity model becomes a two-orbital problem. 
We also assume that the Hund coupling is of Ising type 
since simulating the SU(2) symmetric Hund coupling is 
difficult in QMC. The application of new, more sophis- 
ticated algorithms to this end, such as [321 ] . remains an 
important challenge for the future. 

Besides conventional QMC, we employed the PQMC 
method in the present LDA+DMFT calculation, basi- 
cally following Ref. 28, 0, [3(3] for calculating ground 
state expectation values. With an imaginary time dis- 
cretization At — 0.267cV~ 1 , we take C = 20 time slices 




FIG. 2: (Color online) Spectral function of L1V2O4 at /3 = 
10, 20, 40 eV- 1 (T « 1200, 600 and 300 K); U = 3.6, U' = 2.4, 
and J = 0.6 (all units are in eV). The insets show a closeup 
view immediately above Ef (uj = 0). 



for measurement and "P = 65 time slices before and there- 
after for projection. For the remaining imaginary time 
to (3 — 00, we use the non-interacting Hamiltonian with 
a shifted one-particle potential so that we have n = 1 
electrons/site for the a\ g orbitals and n = 0.25 for the 
orbitals. This shift warranties (approximately) the same 
large-r asymptotic behavior as the interacting Hamilto- 
nian. We performed ~ 3 x 10 8 QMC sweeps and used 
the maximum entropy method for calculating the spec- 
tral function A(u>) and the Fourier transformation of the 
Green function from imaginary time to frequencies, i.e., 
from G(t) to G(iu). 

Results. Let us start with the results of conventional 
QMC at finite T. In Fig. [21 we plot the spectral func- 
tion A(ui), using Coulomb interaction parameters which 
are typical for 3d orbitals, i.e., U = 3.6, U' = 2.4, and 
J = 0.6 eV. The qualitative feature of the result for 
P = 1/T = 10 eV" 1 corresponding to T k 1200 K is 
similar to that of the previous LDA+DMFT calculation 
But for = 40 eV^ 1 (T w 300 K) we note the emer- 
gence of a small structure in the a\ g band just above Ep, 
which is absent for T w 1200 K. At the same time, we see 
no noteworthy temperature dependence for the e g band, 
especially around Ep (see the inset). As far as the finite- 
T QMC is concerned, it is not clear whether the small 
structure in the a\ g band becomes a sharp quasiparticle 
peak at lower temperatures. 

To clarify this point, let us now turn to the PQMC 
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FIG. 4: (Color online) Green function for the ai g (a) and e g 
FIG. 3: (Color online) Same as Fig. H but now at T = ( b ) orbitals obtained by LDA+DMFT(PQMC) as a function 
(PQMC), compared to f3 = 40 eV _1 corresponding to T w 300 of r . compared with the result for (3 = 40 eV 1 (T ~ 300 K). 
K. In agreement with experiment, we can see a sharp peak ^ n tne inset, we magnify the region r > 1. 
slightly above Ef (set to zero) in the ai g band. 



results at T = 0, see Fig. [3] Indeed we can see in the 
PQMC spectrum that the small structure just above Ep 
at T « 300 K becomes a sharp peak, which is consistent 
with the experiment i.e., a peak 4 meV above Ep 
whose width is 10 meV. This a\ g band is lightly doped, 
containing n = 0.98 electrons/site. Unfortunately, the 
exact determination of the renormalization factor (Z) 
from A(u>) or the self energy is difficult because of the 
smallness of the structure and fluctuations from itera- 
tion to iteration in the DMFT cycle. However, the peak 
itself is stable as is the behavior of G(r), plotted in Fig. 
HJ The latter shows a very slow decay for large t in 
PQMC which necessitates the existence of a sharp peak 
at small positive energies in the a\ g band. In contrast, 
for T « 300 K, G(r) vanishes exponentially, see Fig[4] 
(a). 

Discussion. Let us now turn to the physical origin 
of the sharp peak in the a% g band. One possible sce- 
nario which was originally proposed in Rcf. is the 
Kondo effect caused by the hybridization between a\ g 
and orbitals on neighboring sites (note that the on- 
site hybridization is absent). However, it is not trivial 
whether the associated (antiferromagnetic) Kondo cou- 
pling is strong enough to survive a Hund's exchange cou- 
pling as large as J = 0.6. 

To single out the effect of the a\ g -e* g hybridization, wc 
perform an auxiliary LDA+DMFT calculation. To this 
end, we first obtain the a\ g and e g LDA DOS from the 
effective 2-orbital Hamiltonian and then do DMFT cal- 



culations with these DOSes without any hybridization. 
We plot the resulting G(j) and A(uS) of the a\ g band in 
Fig. [5] Clearly, the sharp peak just above Ep survives 
switching off the ai g -e g hybridization. Hence, we can 
conclude that the Kondo scenario due to the hybridiza- 
tion with Cg orbitals [ll[ cannot be the microscopic origin 
for the peak in the a\ g band. Actually, besides contribut- 
ing to the doping of the a\ g band, the e g electrons do 
not play a pronounced role and are only weakly corre- 
lated. Their self energy (not shown) is almost constant 
down to very low frequencies ~ 0.01 eV. The constant 
(ImE ~ — 0.14eV) can be explained by non-interacting 
electrons scattering at disordered spins which is an ap- 
propriate description of the a\ g electrons except for the 
lowest energies. 

Hence, let us turn to the a\ g band itself which is not 
exactly half-filled, but lightly doped with n ~ 0.98 elec- 
trons/site. This suggests that the a\ g band is a lightly- 
hole-doped Mott insulator with a very strongly renormal- 
ized quasiparticle because of the nearness to the doping- 
controlled Mott-Hubbard transition. We can compare 
our results with those of Ref. [33| for the single-band 
Hubbard model on the hypercubic lattice. At n = 0.97 
and very low temperature, these results show a sharp 
peak just above Ep [33|, very similar to our LDA+DMFT 
calculations. An important question for this scenario 
of a doped Mott insulator is whether the strong renor- 
malizations can survive the presence of short-range cor- 
relations beyond DMFT. In this respect, the correlator 
projection method indicates that Z does not vanish for 
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FIG. 5: Spectrum (top) and Green function (bottom) of the 
dig band without ai g -e^ hybridization. The sharp peak above 
Ef survives if the hybridization is switched off. 

the filling-control Mott-Hubbard transition in the two- 
dimensional Hubbard mod el [ 34l| . and also the dynami- 
cal vertex approximations [35I ] shows a strong damping 
of the quasiparticle peak in the vicinity of the Mott- 
Hubbard transition due to antiferromagnetic fluctuations 
beyond-DMFT. However, we believe that such effects are 
less relevant for L1V2O4 because of the frustrated three- 
dimensional lattice and because there is no indication 
that the system is close to a magnetic phase transition. 

In conclusion, realistic LDA+DMFT calculations for 
L1V2O4 show a sharp peak for T — > in agreement with 
photoemission experiments and large renormalizations of 
the effective mass. The physical origin of this peak is the 
lightly doping of the a\ g band which is hence metallic but 
very close to a Mott-Hubbard transition. The HF physics 
is not caused by the hybridization between localized a\ g 
and itinerant e g orbitals. Instead the a\ g orbitals play 
both roles simultaneously whereas the orbitals are 
rather passive and not strongly correlated. 
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